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We report measurements of branching fractions for B — * Kir and B — > irir decays based on a 
data sample of 449 million BB pairs collected at the T(4S) resonance with the Belle detector at 
the KEKB asymmetric-energy e + e~ collider. We also calculate the ratios of partial widths for the 
decays B -> Kit, namely R c = 1.08 ± 0.06 ± 0.08 and R n = 1.08 ± 0.08+° q| ,where the first and 
the second errors are statistical and systematic, respectively. These ratios are sensitive to enhanced 
electroweak penguin contributions from new physics; the new measurements are, however, consistent 
with Standard Model expectations. 
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Tests of the Standard Model (SM) can be performed in 
J5-meson decays to Kit and irir final states, which involve 
various interplays between dominant b — > u tree diagram, 
b — > s, d penguin diagrams and other sub-dominant con- 
tributions. In general, direct comparisons of the mea- 
sured branching fractions with the SM predictions suf- 
fer from large hadronic uncertainties within the current 
theoretical framework. However, many of the uncertain- 
ties cancel out in ratios of branching fractions. Previ- 
ous experimental results [J, 0, 0] for the ratios R c = 
2T(B+ -> K + tt )/T(B+ -> K°tt+) = 1.00 ±0.08 and 
R n = T(B° -> K+Tr-)/2T(B a -*• K°tt°) = 0.82 ± 0.08 
Q deviate from the SM expectations within several ap- 
proaches 0, 0, 0]. For example, Ref. [B[ predicts the 
values R c = 1.15 ± 0.05 and R n = 1.12 ± 0.05, which 
are calculated assuming SU(3) flavor symmetry. If the 
differences between these SM expectations and the mea- 
sured values of R c and R n persist with more data, this 
would imply a large electroweak penguin contribution in 
B —>■ Kit decays [1,0,0]. 

In this letter, we report new measurements of the 
branching fractions for B — > K + tt~ , K + ir°, K°tt°, tt + tt~ 
and 7r + 7r° decays with a data sample five times larger 
than that used in our previous study [l[. Recent Belle 
results for B -► KK, B+ -> K°tt+ and B° -► 7r°7r° 
decays have been reported elsewhere 0, [Io| . The re- 
sults are based on a sample of (449.3 ± 5.7) x 10 6 BB 
pairs collected with the Belle detector at the KEKB e + e~ 



asymmetric-energy (3.5 on 8 GeV) collider ll|. The pro- 
duction rates of B + B~ and B°B° pairs are assumed to 
be equal. The inclusion of the charge-conjugate decay is 
implied, unless explicitly stated otherwise. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of a silicon vertex detector (SVD), 
a 50-layer central drift chamber (CDC), an array of aero- 
gel threshold Cherenkov counters (ACC), a barrel- like ar- 
rangement of time-of-flight scintillation counters, and an 
electromagnetic calorimeter comprised of CsI(Tl) crys- 
tals located inside a superconducting solenoid coil that 
provides a 1.5 T magnetic field. An iron flux-return 
located outside the coil is instrumented to detect JsT£ 
mesons and to identify muons. The detector is described 
in detail elsewhere [12j . Two different inner detector con- 
figurations were used. For the first sample of 152 million 
BB pairs (set I), a 2.0 cm radius beampipe and a three- 
layer silicon vertex detector were used; for the latter 297 
million BB pairs (set II), a 1.5 cm radius beampipe, 
a four-layer silicon detector and a small-cell inner drift 
chamber were used (l3j |. 

Primary charged tracks are required to have a distance 
of closest approach to the interaction point (IP) of less 
than 4 cm in the beam direction (z-axis) and less than 0.1 
cm in the transverse plane. Charged kaons and pions are 
identified using dE/dx information from the CDC and 
Cherenkov light yields in the ACC, which are combined 
to form a K-tt likelihood ratio 1Z(K/tt) — Ck/{£>k +Ar), 
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where Ck (At) is the likelihood that the track is a kaon 
(pion). Charged tracks with TZ(K/n) > 0.6 (<0.4) are 
classified as kaons (pions). Typically, the kaon (pion) 
identification efficiency is 83% (90%), and 6% (12%) of 
selected kaons (pions) are misidentified as pions (kaons). 
Furthermore, we reject charged tracks that are consistent 
with an electron hypothesis. Candidate K° mesons are 
reconstructed as Kg — ► 7r + 7r~ decays with the branching 
fraction taken from Ref. [3] • We pair oppositely-charged 
tracks assuming the pion hypothesis and require the in- 
variant mass of the pair to be within ±18 MeV/c 2 of the 
nominal K s mass. The intersection point of the ir + ir~ 
pair must be displaced from the IP [15J. Pairs of pho- 
tons with invariant masses in the range of 115 MeV/c 2 
< M 77 < 152 MeV/c 2 (±3er) are considered as 7r° candi- 
dates. The photon energy is required to be greater than 
50 MeV in the barrel region, defined as 32° < 7 < 128°, 
and greater than 100 MeV in the end-cap regions, de- 
fined as 17° < 7 < 32° or 128° < 7 < 150°, where 
7 denotes the photon polar angle with respect to the 
direction anti-parallel to the e + beam. 

Candidate B mesons are identified by the "beam- 
energy-constrained" mass, Mbc = \J E*^ m /c 4 — p*g / c 2 , 
and the energy difference, AE = Eg — E^ cam , where 
-^beam ^ s the run-dependent beam energy, and Eg and 
p* B are the reconstructed energy and momentum of the 
B candidates in the center-of-mass (CM) frame, respec- 
tively. Events with M hc > 5.20 GeV/c 2 and \AE\ < 
0.3 GeV are selected for the analysis. 

The dominant background is from e + e~ — > qq (q = 
u, d, s, c) continuum events. We use event topology 
to distinguish the BB events from the jet-like con- 
tinuum background. We combine a set of modified 
Fox- Wolfram moments jl6| into a Fisher discriminant. 
A signal/background likelihood is formed, based on a 
GEANT-based [13 Monte Carlo (MC) simulation, from 
the product of the probability density functions (PDFs) 
for the Fisher discriminant and that for the cosine of the 
polar angle of the B-meson flight direction. Suppression 
of the continuum is achieved by applying a requirement 
on the ratio 1Z = £ s i g / (£ s i s + C q q) , where C s i g {L q q) is the 
signal (continuum) likelihood. Continuum background is 
further su ppr essed through use of the i?-flavor tagging 
algorithm [18j , which provides a discrete variable indicat- 
ing the flavor of the tagging B meson and a continuous 
quality parameter r ranging from (for no flavor-tagging 
information) to 1 (for unambiguous flavor assignment). 
Events with a high value of r are considered well-tagged 
and hence are unlikely to have originated from continuum 
processes. We classify events separately as poorly-tagged 
(r < 0.5) and well-tagged [r > 0.5) in data set I and 
data set II and for each category we determine a con- 
tinuum suppression requirement for 1Z that maximizes 

the value of / <J N%? + N%> . Here, denotes 

the expected signal yields based on MC simulation and 



the average branching fractions of the previous measure- 
ments [H,H, El, and N^ p denotes the expected contin- 
uum yields as estimated from sideband data (Mb c < 5.26 
GeV/c 2 and \AE\ < 0.3 GeV). _ 

Background contributions from Y(45) — > BB events 
are investigated using a large MC sample that includes 
events from b — ► c transitions and charmless B decays. 
After all the selection requirements, no b — ► c background 
is found, while a small contribution from charmless B de- 
cays is present at low AE values for all studied modes. 
Due to K — ir misidentification, large B° — > K + tt~ 
and B + — > K + ir° feed-across backgrounds appear in the 
B° — > ir + ir~ and B + — > 7r + 7r° modes, respectively. 

The signal yields are extracted by performing extended 
unbinned maximum likelihood fits to the (Mb c , AE) dis- 
tributions of the selected candidate events. The likeli- 
hood function for each mode is defined as 

i l,k,j 

where N is the total number of events, i is the event 
identifier, I indicates set I or set II, k distinguishes the 
two r regions and j runs over all components included in 
the fitting function: signal, continuum background, feed- 
across, and charmless B background. The variable Ni^j 
denotes the number of events, and P{ k j — ~Pl,k,j (-^bc 
AE 1 ) are two-dimensional PDFs, which are the same in 
the two r regions for all fit components except for the 
continuum background. 

All the signal PDFs {Pi,k,j=sigaal(M hc , AE)) are pa- 
rameterized by smoothed two-dimensional histograms 
obtained from correctly reconstructed signal MC based 
on the set I and set II detector configurations. Signal MC 
events are generated with the PHOTOS [19] simulation 
package to take into account final-state radiation. Since 
the Mbc signal distribution is dominated by the beam- 
energy spread, we use the signal-peak positions and res- 
olutions obtained from B+ -» L> 7T+ data to refine our 
signal MC (the D° — * K + ir~ir° sub-decay is used for 
modes with a tt° in the final state, while D° — > iV + 7r — 
is used for the other modes). The resolution for the AE 
distribution is calibrated using the invariant mass distri- 
bution of high momentum (phab > 3 GeV/c) D mesons. 
The size of the final-state radiation effects can be assessed 
if we take signal PDFs from MC without PHOTOS and 
use these PDFs to extract the signal yields from the sig- 
nal MC with PHOTOS. The extracted yields decrease by 
5.8% for B° -> K+tt-, 9.4% for B° -> tt+tt" and 3.6% 
for B + — > K + 7T° and B + — * 7r + 7r°, respectively. 

The continuum background PDF is described by a 
product of a linear function for AE and an ARGUS 
function, f(x) = x\fl — x 2 exp [— £(1 — x 2 )], where x = 
AfbccV-E^eam • The overall normalization, AE slope 
and ARGUS parameter £ are free parameters in the fit. 
The background PDFs for charmless B decays are mod- 
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eled by a smoothed two-dimensional histogram, obtained 
from a large MC sample. We also use a smoothed two- 
dimensional histogram to describe the feed-across back- 
ground, since the background events have (Mfc, c , AE) 
shapes similar to the signal, except for a AE peak posi- 
tion shift of ~ 45 MeV. We perform a simultaneous fit 
for B° — ► K + ir~ and B° — ► 7r + 7r~, since these two decay 
modes feed across into each other. The feed-across frac- 
tions are constrained according to the identification effi- 
ciencies and fake rates of kaons and pions. A simultane- 
ous fit is also used for the B + — > K + tt° and B + — > 7r + 7r° 
decay modes. 

When likelihood fits are performed, the yields are al- 
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FIG. 1: 

K+7V-, B° -> tt+tt", B+ -> K+tt°, B + -> tt+tt" and 
B° -> if °7r° candidates. The histo grams show the data, while 
the curves represent the various components from the fit: sig- 
nal (dot-dashed), continuum (dashed), charmless B decays 
(hatched), background from mis- identification (dotted), and 
sum of all components (solid). The Mbc and AE projections 
of the fits are for events that have \AE\ < 0.06 GeV (left) 
and 5.271 GeV/c 2 < M bc < 5.289 GeV/c 2 (right). (A looser 
requirement, -0.14 GeV < AE < 0.06 GeV, is used for the 
modes with a ty° meson in the final state.) 



lowed to float independently for each I (set I or set II) 
and k bin (low or high r region). The Mh c and AE 
projections of the fits are shown in Fig. [TJ while Table U 
summarizes the fit results for each mode. The branch- 
ing fraction of each mode is calculated by dividing the 
total signal yield by the number of BB pairs and by 
the average reconstruction efficiency. The calculation of 
this average efficiency takes into account the differences 
between various I and k bins, and sub-decay branching 
fractions. 

The fitting systematic errors are due to signal PDF 
modeling, charmless B background modeling, and feed- 
across constraints. The first and last of these errors are 
estimated from the fit deviations after varying each pa- 
rameter of the signal PDFs or the yields of the feed-across 
backgrounds by one standard deviation. The effects due 
to fake-rate uncertainties are also included in the system- 
atic error of the feed-across backgrounds. The systematic 
error due to the charmless B background modeling is 
evaluated by requiring that AE > —0.12 GeV, since the 
AE values of the charmless B events are typically smaller 
than —0.12 GeV. The above deviations in the signal yield 
are added in quadrature to obtain the overall systematic 
error due to fitting. 

The MC-data efficiency difference due to the require- 
ment on the likelihood ratio 1Z is investigated with B + — > 
D°n + samples. The systematic error due to the charged- 
track reconstruction efficiency is estimated to be 1% per 
track using partially reconstructed D* events. The sys- 
tematic error due to the 1Z(K/tt) selection is 1.3% for 
pions and 1.5% for kaons, respectively. The Kg recon- 
struction and the systematic error is verified by compar- 
ing the ratio of D + — > Kgir + and D + — » K~tt + it + yields 
with the MC expectations. The 7r° reconstruction effi- 
ciency and the systematic error is verified by comparing 
the ratio of D° -> K+tt- and D° -> X+Tr^vr yields with 
the MC expectations. Possible systematic uncertainties 
due to the description of final-state radiation have been 
studied by comparing the latest theoretical calculations 
with the PHOTOS MC [21] . These uncertainties were 



found to be negligible and thus no systematic error is 
assigned due to PHOTOS. The systematic error due to 
the uncertainty of the total number of BB pairs is 1.3% 
and the error due to signal MC statistics is between 0.4% 
and 0.7%. The final systematic uncertainty is obtained 
by quadratically summing all the contributions, as shown 
in Table HJ 

The ratios of partial widths can be used to extract the 
angle (f>3 and to search for new physics [H, 0, 0] . These 
ratios (listed in Tabic IIIip are obtained from the five 
measurements in Table U and the new measurement of 
B(B+ K°ir+) = (22.8±S:? ± 1.3) x lO" 6 described in 
Ref. 0. The ratio of charged to neutral B meson life- 
time, t b + /t b o = 1.076 ± 0.008 [H, is used to convert the 
branching-fraction ratios into the ratios of partial widths. 
The total errors are reduced because of the cancellation 
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TABLE I: Extracted signal yields, product of efficiencies and 
sub-decay branching ratios (B s ), and calculated branching 
fractions for individual modes. The branching fraction errors 
are statistical and systematic, respectively. 



Mode 


Yield 


Eff.xB s (%) 


£(l(T e ) 


K+TT- 


35851™ 


40.16 


19.9 ±0.4 ±0.8 


7T + 7T~ 




37.98 


5.1 ±0.2 ±0.2 


K+7T° 


U93±% 


26.86 


12.4 ±0.5 ±0.6 


7V+7V 


693tg 


23.63 


6.5 ± 0.4 +° i 


K n° 


379^27 


9.17 


9.2 ± 0.7 ±g:f 



TABLE II: Summary of systematic errors, given in percent. 

o„o 



K+TT- TT + TT" K+n° TT+TT K° TT 



Signal PDF 

Charmless B background 
Feed-across background 
1Z requirement 
Tracking 

7Z(K/n) requirement 
Kg reconstruction 
7T° reconstruction 
# of BB 

Signal MC statistics 



±0.2 

+0.0 
-0.2 
+0.4 
-0.3 

±1.0 

±2.0 
±2.9 
0.0 
0.0 
±1.3 
±0.6 



±0.3 ±0.4 ±0.5 



Total 



±4.0 



+0.6 
-0.0 
+2.2 
-2.1 

±1.0 

±2.0 
±2.8 
0.0 
0.0 
±1.3 
±0.4 

+4.5 



+0.0 

-0.9 
+0.7 
-0.6 

±1.3 

±1.0 
±1.5 
0.0 
±4.0 
±1.3 
±0.4 

+4.8 

-4.9 



+0.0 

-5.0 
+2.5 
-2.4 

±1.4 

±1.0 
±1.3 
0.0 
±4.0 
±1.3 
±0.5 



+ 5.4 

-7.3 



+0.3 

-0.4 
+0.0 
-3.0 

0.0 
±1.5 
0.0 
0.0 
±4.9 
±4.0 
±1.3 
±0.7 

+6.7 

-7.3 



TABLE III: Partial width ratios of B — > Kn and tttt decays. 
The errors are quoted in the same manner as in Table U 



Modes 



Ratio 



2r(K + Tv°)/r(K°TY + ) 

r{K + n-)/2Y{K°n ) 

T(K + ir )/T{K°n + ) 

V{ji + tt-)/T(K + tt-) 

r(7r+7r-)/2r(7r+7r°) 

r(7r+ 7 r°)/r(A' 7r ) 

2r(7r+7r°)/r(A' t V + ) 



± 0.06 ± 
± 0.08 



0.08 



+0.09 
-0.08 



1.08 
1. 

0.94 ± 0.04 ±0.05 
0.26 ± 0.01 ± 0.01 

0.42 ± 0.03 to° i 
0.66 ± 0.07 ± 0.05 

0.57 ± 0.04 



h0.04 
-0.05 



of some common systematic errors. With a factor of five 
times more data than that used for our previous pub- 
lished results the statistical errors on the branching 
fractions for all decay modes are reduced by more than 
a factor of 2.3. The central value of the K q tt° branching 
fraction has decreased from 11.7 x 10 -6 to 9.2 x 10~ 6 
and the K + n~ branching fraction has increased from 
18.5 x 10~ 6 to 19.9 x 10~ 6 , resulting in a change in R n 
from 0.79 ± 0.18 to 1.08 ± 0.12. The obtained value of 
R c = 1.08 ± 0.09 is similar to the previous Belle mea- 
surement (1.09 ± 0.19) but is more precise. The errors 
for R n and R c shown here are the sum in quadrature of 
the statistical and systematic errors. These two ratios 
are now consistent with SM expectations [B|, @, 0, Hi ■ 



In conclusion, we have measured the branching frac- 
tions for B — » Ktt and B — > tttt decays with 449 million 
BB pairs collected at the T(45) resonance with the Belle 
detector. We confirm the expected hierarchy of branch- 
ing fractions : B(K°tt+) > B(K+tt-) > B{K+tt q ) > 
B(K tt°) > B{Tr+ir ) > B{ii + tt-) and find no significant 
deviation from SM expectations in the ratios of partial 
widths. We also find that the ratios R n and R c are both 
in good agreement with SM expectations, in contrast to 
early measurements 0, 0, 0] . 
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